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Abstract. Recently the mechanism was found which allows avoidance of the 

cosmological singularity within the semi-classical formulation of Loop Quantum 
^SJ ' Gravity. Numerical studies show that the presence of self-interaction potential of the 

scalar field allows generation of initial conditions for successful slow-roll inflation. In 
^K^J I this paper qualitative analysis of dynamical system, corresponding to cosmological 

QQ ' equations of Loop Quantum Gravity is performed. The conclusion on singularity 

^^ , avoidance in positively curved cosmological models is confirmed. Two cases are 

considered, the massless (with flat potential) and massive scalar field. Explanation of 
^— V . initial conditions generation for inflation in models with massive scalar field is given. 

^^ ' The bounce is discussed in models with zero spatial curvature and negative potentials. 
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1. Introduction 

Over a long period of time cosmologists believed that quantum cosmology would find 
a way to treat initial as well as final singularity. Loop Quantum Gravity (LQG) is at 
present the main background independent and nonperturbative candidate for quantum 
theory of gravity jU E]. Loop Quantum Cosmology (LQC) is the application of LQG 
to an homogeneous minisuperspace environment. Calculations performed within LQG 
established the possibility to avoid singularity j3]. Semi-classical equations within LQC 
were studied recently and some intriguing conclusions were obtained. First of all, the 
Big Crunch does not appear in models with realistic effective potentials |1] . Moreover, if 
initial conditions are chosen at the bounce so that the scalar field rests at the minimum 
of the effective potential, during the subsequent evolution it climbs the potential and 
hence generates initial conditions for successful slow- roll inflation jHl E] • 

These results were obtained by numerical integration of LQC equations and 
although authors claim that their results are quite generic, application of alternative 
methods to the problem is of interest. Qualitative analysis is quite a powerful method 
and it was successfully applied within the framework of General Relativity (GR) (see 
e.g. IZIIHI) and gauge theories of gravity [Oj. In this paper qualitative analysis of semi- 
classical LQC equations, which are presented in the next section, is performed. The 
origin of bounce is discussed in section 3. Simple model of massless scalar field with 
absolutely fiat potential is studied in section 4. Qualitative behavior of cosmological 
model with massive scalar field is considered in section 5. Conclusions follow in the last 
section. 

2. Semi-classical LQC equations 

Near the classical singularity the picture of continuous spacetime breaks down and 
evolution within LQC is described by difference equations. However there is a region 
where semi-classical differential equations are valid but quantum effects are already 
significant. In this region of small distances the geometrical density a~^ acquires 
quantum corrections. Good approximation to eigenvalues of geometrical density 
operator in LQC for large quantization parameter j is given by JU] 

dj{a) = D{q)a-^, (1) 

where q = {a/a^Y, and a1 = i^-lp is the scale where quantum corrections become 
essential. It can be larger than the planckian scale Ip, since the quantization parameter 
j, which must take half integer values, is arbitrary 01 E]- The quantum correction factor 
D{q) is given by 



Di<l) = [^) g'^' {7 [{<! + ir^' - k - l\''^'] - (2) 

-llg [{q + 1)^/^ - sign(g - l)|g - 1^/^] }^ . (3) 
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Above the scale a^ = 0.357Zp the discrete nature of spacetime becomes unimportant and 
semi-classical equations 13 El IH] are valid 
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where V{ip) is the effective potential, ^/^ = is the scalar field derivative with respect 
to time, if = - is the Hubble parameter, Mp is the planckian mass, the dot denotes 
differentiation with respect to time, and k gives the sign of spatial curvature. According 
to (JH) if y > 0, since the function D is nonnegative, the bounce is possible only for 
positively curved spatial geometry, exactly as within GR, and we assume k = +1 until 
the end of section 5. 

The dynamical system which corresponds to cosmological equations (001) for 
k = +1 can be represented as follows: 
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where| D'^ = dD/dq = Z)(a/a*)~^/(2if). This is four-dimensional system for dynamical 
variables {H,ip,a,ip}. Notice immediately, that the scalar field influence on dynamics 
of the system is due to the unique (last) gradient term in (jH)). When the effective 
potential is completely fiat, V = Vq =const, this term vanishes. Moreover, additional 
simplification is possible here, since the scalar field derivative tp can be expressed from 
(JD) and substituted into (jH)) thus reducing the dimension of the dynamical system to 
two. The system can be completely characterized in terms of {H, a} in that case. 

From the other hand, for many types of effective potentials the gradient of the 
potential can be expressed as a function of the potential itself. Consider massive scalar 

2 2 

field with V = ^^ , where m is a mass, and rewrite the gradient term with the help of 

m 
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This expression can be used in order to reduce the dimension of the dynemical system 
until it does not change sign. This is true for the whole evolution of the dynamical 

X In numerical calculations it is convenient to make interpolated tables for the function D{q) and its 
derivative D' . In classical limit (for large q) the former approaches unity while the latter vanishes. 
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system during bounces and inflation, described e.g. in 0. In that case the system can 
be reduced to three-dimensional one with the following dynamical variables {H^tp^a}. 

3. Origin of bounce in semi-classical LQC 

The possibility of bounce emerges within semi-classical LQC since the phase space differs 
essentially from the case of GR for small scale factors due to quantum corrections. The 
function 

Q(a) ^ (-) ^ (12) 

can change sign in (jTj) so the Hubble parameter derivative may vanish. In the case 
of massless scalar field equality Q = 2 determines the singular point of the dynamical 
system. The right hand side of equation (|T^ (see below) should be used in order to 
determine the singular point position with fixed value of the effective potential§. It 
is the saddle singular point on the a axis in the phase space that ensures the bounce 
for solutions of cosmological equations (see fig. [T]left). One can see that the bounce 
is guaranteed for all phase trajectories crossing the a axis above the saddle Mq, since 
the Hubble parameter changes sign from negative to positive and the transition from 
contraction to expansion takes place. The position of the saddle Mq for massless scalar 
field depends on quantization parameter as 

aAfo^ 1.765 -10-2^5. (13) 

With nonvanishing potential it depends on the value of the effective potential very 
weakly. Only for very large j with large Vq the position of the saddle is shifted a 
little towards smaller values. Expression (fT^ shows that the saddle point comes into 
the region of validity of semi-classical LQC (a > a^) for j > 400. This means that 
trajectories laying below the point Mq appear in the region of semi-classical LQC for 
j > 400. These trajectories do not have bounce since they describe transition from 
expansion to contraction. However, these trajectories originate from and come back to 
the rigion below Oj where purely quantum geometric effects should be responsible for 
the absence of singularity [3]. 

Assume below that the saddle singular point Mq lies below the scale a,, or, 
equivalently, Vq < Mp and j < 400. 

4. Scalar field with fiat potential 

Assume in this section that the effective potential of the scalar field is completely fiat, 
i.e. V = Vo=const. Corresponding dynamical system is given by the following equations 

a = Ha. (15) 

§ We do not consider here all singular points of the system H4IH|I . For instance the usual focus at 
if = ip = 0, which is responsible for oscillations of the scalar field after inflation, is ignored. 



Qualitative Approach to Semi- Classical Loop Quantum Cosm,ology 





Figure 1. Left figure: saddle singular point in the section of the phase space which 
guarantees the bounce. Here j — 30000 just to visualize the saddle. For reasonably 
small quantization parameter j < 400 this singular point lies below the scale of 
quantum geometry a^. Right figure: focus singular point in the phase space for masslss 
scalar field. Phase trajectories are closed curves and the cosmological parameters 
experience oscillations. The quantization parameter is j — 300. Dashed lines surround 
the region where semi-classical LQC is valid, i.e. where |iif | < Mp and a > ai. 



Consider first the case of purely massless scalar field Vq = 0. The phase space 
is represented at fig. [T] (right). The only singular point on the axis a in addition to 
Mo is the focus P. Its position dependence on the value of parameter j, determined 
empirically, is 

ap ^ 0.1944^. (16) 

Therefore, ap > a^ for j > 3.5. This means, oscillating behavior of cosmological 
models with massless scalar field found in jH] is very generic. From the other hand, 
ap > 3.89 when j > 400, i.e. when the saddle of fig. Q enters the scale of semi- 
classical LQC. However, this does not change the qualitative picture of oscillations since 
the asymptotes of the saddle are almost horizontal for small j and almost all phase 
trajectories experience a bounce. 

Now assume the value of the potential is nonvanishing, Vq > 0. Corresponding 
phase portraits are shown at fig El Left figure represents the case of GR. The saddle M 
here originates from the fact that Q = 2 in ()14|) in classical limit and the bounce becomes 
possible for sufficiently large Vq. The thick curve surrounds the region where bounce is 
impossible since the scalar field derivative gets complex and solutions are absent. Phase 
trajectories corresponding to bouncing solutions can cross the a axis only between points 
S and M. At the point S the derivative of the scalar field vanishes. Stability of solutions 
of this type were studied in IJTj. It turns out that they are exponentially unstable. 

The structure of the phase space depends on some critical value of the potential Vc, 
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Figure 2. Phase portraits of cosmological models with the scalar field with flat 
effective potential. Left figure corresponds to the case of GR, Vq — 0.05Mp. Central 
figure represents the phase space with focus and saddle for Vb — 0.006Afp. Right figure 
represents the phase space with Vq = 0.05A/p. Thick curves surround regions where 
the derivative of the scalar field ij) is complex and there are no solutions. For central 
and right figure j = 100. Dashed lines again surround the region where semi-classical 
approach is valid. 



SO if Vo > Vc some singular points disappear. The value Vc is empirically 

Vc = 2.128j^\ (17) 

Central figure describes the phase space for Vq < Vc- There are two singular points 
at the axis a, the focus P and the saddle M. The presence of additional saddle Mg 
below Oj still guarantees the bounce for all solutions. There are two types of possible 
solutions: oscillating solutions, when the scale factor at the bounce amin lies between P 
and A^ (the point where the saddle separatrix bounces), or otherwise quasi-exponential 
bouncing solutions. 

Right figure corresponds to Vq > Vc. Both P and M disappear and the bounce 
below 5* now possible for all solutions. Notice the difference between left and right 
figures. Instability of bouncing solutions within GR is due to the fact that the saddle 
separatrices have de Sitter asymptotes like all solutions. Corresponding separatrices 
do not have de Sitter asymptotic within semi-classical LQC and the class of bouncing 
solutions is much wider. 

It is of interest to consider the dependence of positions of M and P on the values of 
j and Vq. Notice, that although there is dependence on j, the qualitative results remain 
the same for a wide range of quantization parameter variation. Assume first j = 300 
and study only dependence of singular points positions on Vq. The result is shown at 
fig. El 
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Figure 3. Left figure: positions of singular points as functions of the scale factor for 
fixed values of Vq. There are three cases: the function H crosses the axis a in two 
points (yellow and black curves), in one point (blue and red curves) and does not cross 
the a axis (green curve). The effective potential grows from yellow to black curves. 
Right figure: dependence of singular points positions on the value of Vo- Above the 
thick curve -0 is complex and there are no solutions. At the point A the singular points 
disappear and at the point B they appear again. The thick curve marks the scale 
factor values for which the scalar field derivative "ip vanishes. 



Left figure here represents positions of singular points on axis a dependence on the 
scale factor for fixed Vq. The potential has minimal value for yellow curve and grows 
monotonically for red, green, blue and black curves. There are three cases: two points 
M and P are present, only one point is left, no singular points are left. This situation 
is additionally illustrated on the right figure. For purely massless scalar field Vq = 0, 
the point M goes to infinity and we come to the phase portrait represented at fig. ^ 
(right). For small effective potentials both singular points appear in the phase space. 
The position of P depends on Vq rather weakly while the position of M is very sensitive 
to it. Above the thick curve the scalar field derivative becomes complex. Although the 
two singular points come to the phase space for large Vq again, they lie in the region of 
complex ip. 

5. Massive scalar field 



Now we are ready to describe the behavior of cosmological models with massive scalar 
field within semi-classical LQC Although the phase space is four-dimensional, as 
mentioned at the end of section 2, we can use (fTTj) to reduce its dimension to three 
until the oscillations of the scalar field after inflation are not considered. 

Since the bounce is ensured by the presence of the saddle Mq below the quantum 
geometry scale Oj it is convenient to consider initial conditions at the bounce where 
H = 0. From the results of the previous section it follows that the picture of bounce 
is completely determined by the value of the scale factor and effective potential at the 
bounce. When V < Vc, where Vc is given by (|T7j) there are two singular points in the 
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section of the phase space {H^a} and oscillations develop if the scale factor lies between 
the points P and A^. However, during oscillations the value of the scalar field increases 
(if '0 > 0) or decreases (when ■?/' < 0) and the scalar field climbs the potential walls 
jnj. So in the section {H^a} of the phase space the singular points P and M tend to 
each other until the value of the potential reaches the critical value Vc- When V = Vc 
possibility of oscillations disappear and the phase trajectory inevitably escapes from the 
bounded region near the point P. The scale factor starts to grow quasi-exponentially, 
and the solution approaches classical region so inflationary regime of GR develops. When 
inflationary expansion begins the kinetic energy of the scalar field becomes suppressed 
and the subsequent expansion is determined by slow-roll of the scalar field. Notice that 
the less is the scalar field mass (the more fiat is the potential), the more oscillations can 
develop before the infiationary stage begins. When V > Vc there are no oscillations and 
the bounce is unique. 

The same qualitative behavior should take place for all symmetric nonnegative 
effective potentials of the scalar field. It is interesting to note, that the possibility of 
bounce for zero curvature models {k = 0) appears in semi-classical LQC with negative 
potentials. The phase portrait of zero curvature model with Vq < is similar to fig. ^ 
(right). 

Notice at last, that the validity of the phase analisys presented here is limited. For 
V > Mp for instance the quantum nature of the scalar field becomes important. Also 
the Hubble parameter is constrained by the usual condition H < Mp. Finally, the 
requirement a > ai should be satisfied. 

6. Conclusions 

Qualitative study of dynamical system corresponding to semi-classical LQC equations 
is performed. It is shown, that the bouncing character of cosmological solutions within 
semi-classical LQC is guaranteed by the presence of the saddle singular point on the 
scale factor axis. The bounce takes place for all solutions already in semi-classical 
region of LQC if j < 400. Although some 'singular' solutions appear in semi-classical 
region for j > 400 they should be free of singularity due to quantum descrete structure 
of spacetime [3] which is not considered here. The possibility of oscillations before 
the bounce in models with massive scalar field is explained. It is shown that initially 
oscillating solutions inevitably followed by the infiationary stage. The bounce is possible 
also for zero curvature models within semi-classical LQC with negative potentials. 
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